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Abstract 

We show that for 1-1-1 dimensional Causal Dynamical Triangulations (CDT) 
coupled to 4 massive scalar helds one can construct an effective transfer matrix if 
the masses squared is larger than or equal to 0.05. The properties of this transfer 
matrix can explain why CDT coupled to matter can behave completely different 
from “pure” CDT. We identify the important critical exponent in the effective 
action, which may determine the universality class of the model. 
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1 Introduction 


Causal Dynamical Triangulations (CDT) in 1+1 dimensions can be considered a 
toy model for more advanced models of quantum gravity. The simplest version of 
the model, pure fluctuating geometry without matter, can be solved analyticalljj^ 
[2]. The random geometries present in the CDT path integral are constructed 
by gluing together flat simplices (triangles) in such a way that one has a global 
time foliation. The topology of space at a given time is assumed to be S^. This 
topology is preserved in the time evolution. In all our considerations time is Wick 
rotated and the triangles used are assumed to be equilateral with edge lengths a. 
We label time with integer numbers t. The geometry of a spatial slice at time t is 
completely characterized by its length, i.e. in the CDT model by n{t) - the number 
of edges forming the spatial . In the path integral representation of the time 
evolution, spatial states at integer times t and t + 1 are connected in all possible 
ways consistent with the foliation. In the case of CDT without matter fields the 
time evolution is generated by a transfer matrix {nt\M.\nt+i) = exp(—L(nt, nt+\)) 
with correctly normalized spatial states |n)[^ The explicit expression for L{n,m) 
can be interpreted as a term of the effective action, since in obtaining it we sum 
over all geometrical realizations joining the two states. 

One may view the same geometry using a dual trivalent lattice, where vertices 
are located in the centers of triangles and links being dual to the edges in the 
original lattice. Each vertex has exactly three neighbors, two at the same time 
t (which can be considered as a half-integer time with respect to the original 
triangulation) and one at time f ± 1. As before the links at the same time value 
form a closed spatial geometry S^, the quantities rit and ut+i represent numbers 
of links pointing up or down from the line at f + 1/2. The dual formulation is 
completely equivalent to the original one and we will use it in this article. 

The Hilbert-Einstein action for a triangulation T, S{T), provides the weight 
exp(—S'(T)) to be assigned to the triangulation in the path integral. In 1+1 
dimensions there is no curvature term (it is a topological invariant), and we are 
left with the cosmological term, which on the lattice takes to form 

S'(T) = ^ Xrii = XN, rii = (1) 

i 

where N is the total number of triangles (dual vertices) of the triangulation T. 
A is the dimensionless bare cosmological constant. The weight contains a factor 

e-¥^+m) = gn+m^ ^ = 6"+ (2) 

^There are a number of generalizations, also only dealing with fluctuating geometries, which 
can be solved analytically [T]. 

^The normalization should include the symmetry factor of the states, in the original paper 
they were realized as marking one of the triangles in {n{t)\. It is possible to take a symmetric 
definition, which we do in this paper. In the following we assume the norm {n\m) = Sum- 
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in each transfer matrix element (n|A^|m). This factor compensates the entropy 
factor, always present in the models of quantum geometry, where the number of 
different triangulations between n and m typically (for large n and m) grows as 
exp(AciV) with some critical Ac- The quantum amplitude (the partition function) 
can now be written as 

Z = ^exp(-^(r)) = exp ^-^L(nt,nt+i)j (3) 

In [3]-[5]) we studied the effect of coupling d massive scalar helds to the 1+1 
dimensional CDT model. We assumed the helds to be located at the vertices of 
the dual lattice and introduced the partition function 

Z = J2 [ll^€ - Smatterim) ( 4 ) 

r '' 

with the Gaussian matter action 

Smatter = , ( 5 ) 

i^fi 

where Uj denote the links in the dual lattice and p = 1,... ,d. 

It is clear that the integral over the helds (|^ is well dehned for > 0. In 
order to see the role of the mass parameter in the action Smatter, we redehne the 
held variables by = mcf)^, and the matter action changes to 

Scatter = ^ ~ ^ 

Ij^j ^fi i,fj, 

From (|^ it follows that the mass controls the range of held interactions: in the 
large mass limit —)■ oo, the couplings between neighboring vertices can be 
neglected, and as a consequence, the contribution of matter can be eliminated 
{d = 0, pure gravity), while in the small mass limit —)■ 0 (the massless case), 

the range of interaction becomes long. It is also clear that there may no longer 
exist a transfer matrix depending only on the geometric variables n* at time slice 
t and its neighboring time slices f ± 1. Integrating out the matter degrees of 
freedom might introduce long range interactions between various time slices, an 
ehect which one would expect to increase with decreasing the mass. Indeed, the 
ehect of small or zero mass matter helds on the 1+1 dimensional global geometry 
is dramatic when d > 1 as reported in [3] (and also seen in earlier studies using 
other kind of matter helds [6]). One observed the appearance of a “semi-classical” 
de Sitter-like blob with Hausdorh dimension dn = 3, much like what has been 
observed in higher dimensional CDT [7] (see [8] for a review). Surprisingly, even 
for this system there seems to exist an effective local action which couples only 
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neighboring time slices and which describes the “semi-classical” blob and the 
fluctuations around it. 

Inspired by these results we will try clarify to which extent we for d > 0 can 
talk about an effective transfer matrix depending only on the geometric variables 
Tit, and try to understand which characteristic features of the transfer matrix 
change when a semiclassical blob is created for d > 1. 

More precisely we will study 1-1-1 dimensional CDT with d = 4 and = 0.05. 
This mass is so small that a blob is formed and still so large that an effective 
transfer matrix depending only on rit can relatively easily be extracted (for smaller 
masses this becomes increasingly difficult). We will compare our results to that 
of 1-1-1 dimensional CDT without matter fields where no blob is formed. 

The form of the transfer matrix for pure 1-1-1 dimensional CDT will be an 
indication of the form to be expected for the other cases. It is laiH] 


{n\M\m) 



n + mY 


^n-\-m _ -L(n,m) 

y ^ 


(7) 


We are interested in the asymptotic limit, when n, m are large but (n — m) stays 
finite. We obtain 


1 fu — iriY 

L{n,m) = C — {n + m) \og{2g) + - - \-V{n + m) 


2 n -|- m 


where the potential 


V(n + m) = ]- log(n + m) + —— - -r 0(- 

^ ' 2 ^ ' A{n + m) b 


n + mY 


) 


( 8 ) 

(9) 


In ([^ we see a term —(n + m) log( 25 f) = —(n + m) log{g/gc), a linear term related 
to the entropy of states. A similar entropy term will always be present and we 
will fine tune the value of g to be as close to gc as possible in numerical analysis 
aimed at the determination of the transfer matrix elements. The next term is 
a “kinetic” term, coupling the spatial volumes at slices t and t -|- 1. This term 
is “local” and we also expect such a term to be present in the effective transfer 
matrix. Finally we have a “potential”, diagonal term V(n + m). The leading 
large volume term in D is a logarithmic term. The rest of the terms decrease for 
large spatial volume. One can show that if we observe a blob with dn = ^ one 
cannot have terms of the type (u m)" with 0 < a < 1 in D. As a consequence 
the value of the parameter in front of the potential log(n -|- m) is important for 
the global behavior of the model, as we will explicitly show later. The small- 
volume corrections may be important for a detailed behavior of the system at 
small volumes. We will not be concerned with the detailed analytic form of these 
corrections for CDT coupled to matter, but will determine them numerically by 
Monte Carlo simulations. More precisely we parameterize the transfer matrix 


{n\M^^\m) = exp(—Le//(n, m)) 


( 10 ) 
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as 


Leff{n, m) = C - {n + m) \og{g/gc) + /i log(n + m) + 


1 {n 


m] 


r (n + m) ’ 


( 11 ) 


for (n + m) > K for some K, while for [n + m) < K the transfer matrix is simply 
determined from the compnter simulations. In an overlap region we match the 
assumed, parameterized transfer matrix to the numerically determined one and in 
this way we determine the coefficients /x and F, which will depend on the number 
of Gaussian fields, d, and the value of the mass squared, m?. Once the constants 


in Lf,ff are determined we can use the transfer matrix (10) for arbitrary large 
n + m > K. 

If a transfer matrix {n\Ai\m) exists also for systems with matter coupled to 
the geometry then the partition function for a system with periodic boundary 
conditions in the time variable with period T is given by 


Z{g,T) = '^{ni\M\n 2 ){n 2 \M\n‘i) ■ ■ ■ {nT\M\ni) = ti, ni = n{ti) (12) 

{rii} 

This function is defined for g < g^ and the limit g ^ gc corresponds to taking 
a large volume limit. Essentially, Z{g,T) is the partition function we will use 
when checking that the transfer matrix model produces the same spatial volume 
distributions as the full model defined by (|^. 

The rest of this article is organized as follows. In the next Section we discuss 
the numerical methods used to obtain the estimate for gc and to determine the 
elements of the transfer matrix. To see if the transfer matrix determined this 
way reproduces the observed distributions of spatial volumes we perform Monte 
Carlo volume experiments, using the determined transfer matrix for systems with 
a fixed volume N and periodicity T, comparable to that used in the paper [3] . In 
Section 3 we analyze the eigenvalue spectrum of the transfer matrix for g ^ gc 
for pure CDT and for CDT with 4 scalar massive field with mass rn? = 0.05. We 
recall the exact dependence of the spectrum for pure gravity and compare it with 
the d = 4 case. We summarize the results in section 4. 


2 Determination of the transfer matrix 

As described above we assume that after integrating out the matter fields in (Q 
that the quantum geometry can be described by an effective transfer matrix Ai 
with matrix elements {n\Ai\m). The matrix elements are (semi-)positive and the 
matrix is symmetric. The last property follows from the invariance of the model 
under the change of time arrow. We will try to determine the effective transfer 
matrix from Monte Carlo simulations. We use the numerical set-up described in 
[311115]. For a periodic system with the period T the probability to observe a 
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sequence of spatial volumes {ni, n 2 ,... ut} is 


{ni\M\n2){n2\M\n^) ■ ■ ■ {nT\M\ni) 

"2.... nr) =-- 

This quantity is invariant under cyclic permutations and reversion of order 
ni, 77-2,..., riT —)■ riT,..., 77 , 2 , ni. The probability to observe a particular combina¬ 
tion of spatial volumes can be measured by Monte Carlo simulations. Using such 
measured probabilities for a sequence of periods T we can determine {n\Ai\m) 
up to a multiplicative constant. We compare two measured probabilities 



PT,i{n,m) 


(n IW11 777 ) (777 1W1 1 77 ) 
tr 


(14) 


where we determine the probability to observe volumes n and m in two neigh¬ 
boring slices for a system with the period T -|- 1 and 


P2T,T{n,m) 


( 77 1^ 1777 ) ( 777 1^ 1 77 ) 

tr 7V12T 


(15) 


where we determine the probability to observe volumes n and m in two slices 
separated by T steps for a system with the period 2T . We determine the transfer 
matrix elements from 


(77|A1|777) = C 


PT,i{n,m) 

^/P2T,T{n,m) 


(16) 


where C is the multiplicative factor mentioned above. Notice that the Ihs of (16) 
is independent of T. 

The transfer matrix depends on the parameter g and we expect the behavior 
(77|W1|777) oc ((y'/s'c)"'"''™'- We are interested in the for g —>■ g^, where the average 
volume diverges (the continuum limit). To avoid the problem of infinite volume 
we introduce in the simulations a modified action 


T 

S{T) ^ S{T) + e - ^o)^- 

7=1 


(17) 


This modihcation forces the spatial volume at each slice to fluctuate around uq. 
The new action leads to a modihed transfer matrix but one can easily recover the 
original transfer matrix, as will be explained below. The modification of the action 
serves two purposes. Firstly, it allows us to determine the critical coupling gc. The 
average volume distribution in the slices, for T and no sufficiently large, should 
have a maximum at { 77 ) = no- For finite T we observe different distributions for 
different choices of no, but for g = gc the maximum should be stable. Typical 
distributions are shown in Fig. for a system with four massive scalar fields 
with a mass 777 ^ = 0.05 and g = gc- This method can be used to determine the 
critical parameter gc even in cases, where the transfer matrix may not be a good 
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Figure 1: Averaged volume distribution in a slice for a system with d = 4 and 
= 0.05. We take Uq = 100, 200 and check the location of a maximum using 
T = 8 and T = 14. 
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m 

Qc estimates 

the best estimate Qc 

0.00 

0.2988-0.3000 

0.2991 

0.05 

0.3205-0.3212 

0.3210 

0.10 

0.33375-0.3338 

0.33378 

0.15 

0.34412-0.34417 

0.34415 

0.20 

0.35284-0.35285 

0.35285 

5.00 

0.48322-.048342 

0.48341 

Pure gravity 

0.4999-0.50 

0.50 


Table 1: Estimated values of Qc- The estimates come from volume distributions 
with no = 0, 50, 100, 150, 200 and T = 20. 

approximation. In the table [T] we list the estimates for a range of mass parameters 
m? and, for comparison, for pure gravity. In the following we use g = Qc- 

Secondly, the modihcation ( [ItI ) makes the spatial volume fluctuations much 
more controllable in the Monte Carlo simulations. It changes the probability 
distributions to 

~ {ni\M\n2){n2\M\n^)---{nT\M\ni) 

Tr(ni,n2,...,nT) =-TATT- 

where the relation between Ad and AA is given by 

(n|M|m) = (19) 

again up to an arbitrary factor C. This relation permits us to determine the 
matrix elements {n\Ai\m). The method will work for n, m restricted to a window 
around no, where volume fluctuations are small. The size of the window depends 
on the parameter e and typically outside this window the fluctuations become 
large. In practice we measure the matrix elements for a sequence of no and 
connect the results by requiring the best overlap. 
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2 

m 

0.05 

0.10 

0.15 

0.20 

5.0 

Pure Gravity 

/i 

0.64628 

0.59365 

0.55432 

0.53233 

0.503512 

0.505763 

r 

2.5 ±0.05 

2.2 ±0.1 

2.17 ±0.15 

2.15 ±0.1 

2 ±0.1 

2 ±0.05 


Table 2: Estimated values of the parameters /i and T as explained in the text. 


The transfer matrix, (lO)-(ll), depends on the two parameters /i and T. We 
determine /i by measuring the diagonal part of the transfer matrix (i.e. {n\M.\n)), 
where the kinetic term vanishes. We use T = 8. For T > 5 the results seem 
insensitive to the periodic boundary conditions imposed in the time direction. 
As before we combine the measurements for various overlapping windows around 
different values of uq. The results of best hts are shown in Fig. for pure gravity 
(as a test of the method) and for d = 4, = 0.05. 


1.00 
0.95 
^0.90 
§0.85 
V 0.80 

ao 

o 

0.75 

0.70 

0.65 




Figure 2: The logarithm of the diagonal of the transfer matrix aX g = Qc for pure 
gravity (left) and for d = 4, m? = 0.05 (right). The plots show a combination of 
results obtained with uq = 0, 50, 100, 150. The black line represents a £t to the 
asymptotic power behavior of the diagonal part, with the coefficient determined 
in the range limit by the blue dots. 


We determine T by measuring the transfer matrix for n + m hxed (such that 
only the kinetic term changes). For sufficiently large c we have 


(n|M|c — n) = N‘{c) exp 


(2n — cY 
fc 


( 20 ) 


where the terms in the effective action which only depend on c are included in 
the normalization. This is illustrated on Fig. and we observe that the width 
of the Gaussian grows with c. This growth is with a very high accuracy linear, 
as is illustrated on the hgure The fitted values of the parameters /i and T are 
listed in the table |2l 

As a check of the quality of the transfer matrix determined this way we per¬ 
formed a Monte Carlo simulation where the probability assigned to a geometry 






















Figure 3: The transfer matrix {n\M\c — n) for various c's: pure gravity (left) and 
d = 4, rn? = 0.05 (right). 



Figure 4: Linear behavior of k{c) ~ Fc for pure gravity (left) and d = 4, m? = 0.05 
(right). 


with spatial volumes {ui, n 2 ,..., ut} is 
F(ni, 77.2,..., Ut) oc ^\n2){n2\M‘^^^\n^) ■ ■ ■ (nyl 

( 21 ) 


The “effective” transfer matrix entries {ni\Ai^^^\n 2 ) used in (21) are 


m) = 


(771 1777.), m + n < K 

(77 1 Al*^ 1 777), m + n>K. 


( 22 ) 


The small volume part is obtained numerically and the large volume part is given 
by (10) and ( [II| ) with /7 and F determined as described above. K is chosen in 
the range 50-100. The extra factor in (21) is added to enforce the total volume 


to fluctuate around a given value N. The resulting distribution can be compared 
to the one obtained using the full partition function (|^. We use here the same 
method to analyze spatial volume profiles as that presented in 13 a E]: each 
configuration is shifted in such a way that the “center of mass” of the spatial 
volume distribution is placed at time t = T/2. In this way we produce an artificial 
maximum also for pure gravity where there is no blob, but the properties of this 
distribution is quite different from the “real” blob distribution, as explained in 
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EISIE] and as seen from Fig. The agreement is very good. For the case 




Figure 5: Comparison of n{t)/N using the effective transfer matrix and full CDT: 
pure gravity (left) and d = A, m? = 0.05 (right). We use N = 16000, T=200 for 
both cases. 


d = 4, = 0.05 the small volume part of the transfer matrix is important 

in order to get a quantitative agreement of the volume prohle in the tail of the 
distribution (the agreement becomes less good at a quantitative level if iF < 50 
in (22)) The existence or non-existence of the blob is however entirely linked to 
the value of jj, in the effective action (0: simply using this effective action (with 
g = gc) and adding a small volume term l/(n -f m) like in (|^ in order to stabilize 
the logarithmic term for n -|- m = 0 one obtains by Monte Carlo simulations a 
transition from a non-blob phase to a blob phase simply by changing fi. /i is like 
a critical exponent: it is well know that /i = 1/2 appears as an entropy in CDT 
[2], and this is seen explicitly from ([^ for large n\ 

{n\M\n) ~ (2n)-‘/"(g/g,)^”(l + 0(l/n)), (23) 


Changing the details of the triangulations used in CDT, changing the scale of n 
etc, will change gc but not the exponent /x. Similarly, when we couple CDT to 
the matter helds we have 


{n\M\n) {2n) ^^{g/gcf'^il + 0{l/n)). (24) 

where /x will invariant under changes of triangulation details, but will depend on 
d and m. Our numerical results suggest that there is a critical value /Xc such that 
for /i > /Xc the geometry change universality class from that of pure CDT (with 
Hausdorff dimension dn = 2) to that the “blob” geometry (which has djx = 3). 


3 Eigenvalue spectrum of the transfer matrix 


The analytic solution of H-l dimensional CDT |2] permits us to determine the 
eigenvalue spectrum of the exact transfer matrix as a function of g. Using the 
parametrization 

1 _ 1 

^ “ 2 cosh(/3) ’ 
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and solving the eigenvalne eqnation we find (see also [0] for the eigenvalne spec- 
trnm for a more general model) 

A„ = n = l,2 ,... (26) 

The effective transfer matrix determined above for d = 4 was only obtained np 
to a normalization so we can only determine the ratio of eigenvalnes for this “em¬ 
pirical” transfer matrix. For the pnre CDT model we have A„/Ai = exp(—2?7,/3). 
Thns the ratio goes to zero exponentially with n for g < Qc- For g = Qc the 
spectrnm becomes degenerate. In a nnmerical analysis we will never be able to 
achieve this limit, since onr nnmerical transfer matrix is necessarily finite and the 
dependence on Umax for g = gc for pnre gravity is illnstrated on the figure We 
can see that at the critical point the convergence of the ratios to one is slow as 
function of rimax- 



Figure 6: Left: Eigenvalues for pure gravity at g = gc with cut-off rimax = 
200,400,600,800. Right: A 2 /A 1 as a function of Umax ■ 

We repeat the same analysis for the case d = 4, = 0.05. The dependence 

of the eigenvalue spectrum at the critical point for various cut-off values rimax is 
shown in Fig. We see that even at the critical point the first few eigenvalues 
become cut-off independent. This can be attributed to a faster fall-off of the 
transfer matrix elements for large volumes. 

The eigenvalue spectrum is markedly different from the pure CDT case since 
the there is a gap between the first and the second eigenvalue even at g = g^ a gap 
with does not vanish for large Umax- On the other hand the rest of the eigenvalues 
behave like the pure CDT eigenvalues in the sense that they coincide for growing 
^max- The separation of the largest eigenvalue from the rest is a reflection of the 
existence of a stalk. In fact let us denote the first eigenvector (with the largest 
eigenvalue) ni(n). Fig. shows that for small n we have z/i(n)^ = P{n), where 
P{n) is the probability distribution for spatial volumes in the stalk (which is 
almost independent of T and N). 
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Figure 7: Eigenvalues for d = 4 m? = 0.05 and cut-off Umax = 200, 300, 400, 600. 
The first four eigenvalues are essentially independent of Umax- 



Figure 8: The distribution of spatial volume in the “stalk” for <7 = 4 and m? = 
0.05 and the square of the hrst eigenvector of the transfer matrix. 

4 Conclusions 

In CDT there exists a transfer matrix. This transfer matrix depends on the 
geometry of the spatial slices. In 1-1-1 dimensions the geometry of the spatial slice 
is fully characterized by its length (assuming the topology is that of S^). Once 
we couple the geometry to matter the transfer matrix still exists, but it will also 
depend on the matter degrees of freedom. Integrating out the matter degrees of 
freedom might introduce non-local interactions and invalidate any simple transfer 
matrix description in terms of geometry only. However, it turned out that for 
massive free Gaussian helds coupled to geometries, and for the mass not too small 
there is such an effective transfer matrix which describes very well the fluctuating 
geometry of the full model. 

We determined the effective transfer matrix in 1-1-1 dimensional CDT coupled 
to 4 massive scalar helds with > 0.05. 0.05 was the smallest value of 
where we could reliable determine an effective transfer matrix. We found that 
the most important term in the effective transfer matrix was an “entropic” factor 
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l/(n + m)^. /X is like a critical exponent, much like the entropy or susceptibility 
exponent 7 in non-critical string theory or the theory of dynamical triangula¬ 
tions. In the case of non-critical string theory 7 depends on the matter coupled 
to the 2D geometry and there is a phase transition between two completely dif¬ 
ferent classes of geometries at 7 = 0. For 7 > 0 the two-dimensional geometry 
degenerates into so-called branched polymers. We have a somewhat similar sce¬ 
nario here: /i depends on d and m? and there exists a /Xc such that for ^ 
the geometry undergoes a phase transition and develops a “blob” with Hausdorff 
dimension du = 3. The appearance of the blob had a profound impact on the 
effective transfer matrix. A gap developed between the two largest eigenvalues 
of the effective transfer matrix and the eigenvector corresponding to the largest 
eigenvalue was essentially equal to the square root of the probability distribution 
of spatial volumes of the stalk associated with the blob. We conjecture that a 
similar effective description of the blob-non-blob dynamics will be present for 
higher dimensional CDT where it has been shown that there also is an effective 
transfer matrix which describes distribution and fluctuation of the spatial volume 
of the time slices HD]. 
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